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ABSTRACT 

An analysis of the kinetics of H 2 formation on interstellar dust grains is presented 
using rate equations. It is shown that semi-empirical expressions that appeared in 
the literature represent two different physical regimes. In particular, it is shown that 
the expression given by Hollenbach, Werner & Salpeter [ApJ, 163, 165 (1971)] applies 
when high flux, or high mobility, of H atoms on the surface of a grain, makes it 
very unlikely that H atoms evaporate before they meet each other and recombine. 
The expression of Pirronello et al. [ApJ, 483, L131 (1997)] - deduced on the basis of 
accurate measurements on realistic dust analogue - applies to the opposite regime (low 
coverage and low mobility) . The implications of this analysis for the understanding of 
the processes dominating in the Interstellar Medium are discussed. 

Key words: ISM:abundances - ISM: molecules - ISM:atoms - atomic processes - 
molecular processes 



1 INTRODUCTION 

The problem of the formation of molecular hydrogen, the 
most important species in the universe, is a fundam ental 



chowski ( |1981|; 1983| ); Aronowitz & Chang (|L985|); Duley & 
Williams (^986; ); Buch & Zhang ( |l99l| ); Sandford & Alla- 
mandola ( |1993| )] than experiments have been done. 

A milestone has been certainly set by Hollenbach & 



open question in astroph ysics ( Duk^fc^ Williams 1984 ). It 
was recognised long ago (Glould & Salpeter 1963) that H2 
cannot form in the gas phase in the Interstellar Medium 
(ISM) efficiently enough to account for its abundance. It was 
proposed that dust grains act as catalysts allowing the pro- 
tomolecule to quickly release the 4.5 eV of excess energy (in a 
time comparable to the vibration period of the highly vibra- 
tionally excited state in which it is formed). Briefly stated, 
the problem is as follows. An H atom approaching the sur- 
face of a grain has a probability £ (sticking coefficient) to be- 
come trapped. The adsorbed H atom (adatom) will spend an 
average time tn (residence time) before leaving the surface. 
If during the residence time the H adatom encounters an- 
other H adatom (which has just landed on the surface or was 
already trapped in a deeper adsorption site), an H2 molecule 
will form with a certain probability. Given the fact that, un- 
til recently, there were no experiments done in conditions 
relevant to the ISM, and that little is known of the chemical 
composition and morphology of dust grains, it is not surpris- 
ing that quite different models co-existed. This is an area, 
as many others in astrophysics and astrochemistry, in which 
by far more theo retical pap ers ha ve been written [Gould 
fc Salpeter ( |l963| ); Williams ( ]l968[ ) ; Hollenba ch fc Salpeter 
(|l970|; |l97l|); Hollenbach, Werner & Salpeter (|l97l|); Smolu- 



Salpeter (1970; 1971), who treated sticking and accommo- 
dation of H atoms in a semiclassical way, while the mobil- 
ity was treated quantum mechanically. They concluded that 
tunneling between adsorption sites, even at 10K, would have 
assured the required mobility. Hollenbach et al. (1971) ob- 



tained for the steady state production rate of molecular hy- 
drogen per unit volume the simple expression: 



Rh 2 



1 



riHVHO-^rjrig, 



(1) 



where nn and «h are the number density and the speed 
of H atoms in the gas phase respectively, a the average 
cross-sectional area of a grain, n g is the number density of 
dust grains, and rj is the probability that two H adatoms 
on the surface meet and recombine to form H2. Not e tha t 
in the original formulation by Hollenbach et al. (1971) 
Rh 2 — (l/2)nuVHC T "/n g , where 7 is the fraction of H atoms 
striking the grain that eventually form a molecule, namely, 
7 = £77. Eq. (nf) states that, for r} = 1, whenever two H atoms 
are adsorbed on a grain, a H2 molecule is formed. 

From the experimental poi nt of view, besides some pio- 
neering work in the e arly '60s (King & Wise 1963) and '70s 
( [gchuttc ct al. 197tj ), o nly very recently the problem has 
been investigated again ( Pirronello et al. 1997a ; Pirronello et 
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al. 1997b). Pirronello et al. performed their measurements in 
an ultra high-vacuum (UHV) chamber (typical experimen- 
tal pressures were in the 10 -10 torr) irradiating the sample, 
maintained at temperatures between 5K and 15K, with H 
and D atoms from two different triple differentially pumped 
lines (D atoms were used to obtain a better signal to noise 
ratio). For the very first time they used as a substrate a nat- 
ural "olivine" (a Mg, Fe silicate) slab (mechanically polished 
until shiny), that has to be considered with good reasons a 
better analogue of interstellar dust than any model surface. 
A quadrupole mass spectrometer detected the amount of HD 
formed on the cold substrate. Measurements were performed 
both during and after irradiation with H and D atoms. In the 
latter case, a Temperature Programmed Desorption (TPD) 
experiment was carried out in which the temperature of the 
sample is quickly ramped to over 30 K to desorb all weakly 
adsorbed species. 



The main results obtained by Pirronello et al. (1997a 



1997b) are: (a) In the temperature range of interest for inter- 



stellar applications (between 10K and 15K), the formation 
rates deduced from their experimental data are up to one 
order of magnitu de lower th an those calculated by Hollen - 
bach & Salpeter fll97C| ; |l97l[ ) and Hollenbach et al. ( |l97l| ); 

(b) According to their desorption spectra, hydrogen, that 
is adsorbed as atomic, becomes mobile only around 10 K, 
even at the h igh coverage regime [see Fig. 2 in Pirronello 
et al. ( 1997b| )]. Thus, at the lowest temperatures (less than 
about 10 K) tunneling alone does not provide enough mobil- 
ity to H adatoms. Instead, they find that thermal activation 
is required. A possible scenario is that thermal energy is 
necessary to raise H adatoms inside the adsorption well to 
an energy level from which tunneling can become effective; 

(c) According to a careful analysis of the kinetics of HD 
desorption spectra during TPDs, different regimes can be 
recognised during H2 formation, depending on the values of 
certain parameters discussed below. This analysis has led to 
the proposal of the following expression for the steady state 
formation rate of H2, 



-Rh 2 = -(nH«Htf£t H ) n s N vf(T,a,6E)jt, 



(2) 



where JV 2 is the average number of hops an adatom needs 
to make to encounter another adatom while performing a 
random walk, and vf(T, a, SE) describes the hopping rate 
of adatoms due to both thermal activation and tunneling. 
Here, v is a characteristic attempt rate, T is the grain tem- 
perature, while a and SE are width and height of the en- 
ergy barrier, respectively. 7/ is the probability that two H 
adatoms recombine after encountering. In this expression 
the rate of H2 formation is proportional to the square of the 
effective incoming H flux (nH«HO"C) on the grain surface, and 
is built on a purely phenomenological basis to interpret the 
experimental data on desorption kinetics. 

_In Eq. (Flj), that has been used by Hollenbach et al. 



(1971) and others in their chemical models, the recombi- 



nation rate is linearly proportional to the effective incoming 
H flux. In this note it is shown that rate equations yield Eqs. 
(jjj) and as two distinct limiting cases of the H recombi- 
nation rate. 



2 THE CALCULATIONS 

For simplicity we will consider the H2 production rate on a 
single grain; the total production rate per unit volume can 
be obtained by multiplying it by n g , the number density of 
dust grains. The number of H adatoms on the surface Ni, 
and that of H2 molecules N 2 , as a function of time is given 
by the rate equations 



Ni = F - piNi - 2a.Ni 
N 2 = aN 1 2 - P2 N 2 



(3a) 
(3b) 



where F is the rate of adsorption of H atoms on the grain, a 
is the rate of H2 recombination (given by the product of the 
diffusion coefficient and the recombination probability), pi 
and p 2 are the desorption rates for H and H2, respectively. 
The first term in Eq. ( pa| ) represents the growth of the H 
population due to the incoming flux; the second term is the 
decrease in the H population due to the desorption of H 
adatoms; and the third term represents the rate at which 
H adatoms are lost due to the recombination process. In 
Eq. (3b) the first term represents the rate of creation of H2 
molecules. It is related to the last term in Eq. (^) through a 
factor 1/2 because two H adatoms are needed to form one H2 
molecule. The second term in Eq. (3b) is the desorption rate 
of H2 molecules and is equal to the desired production rate 
Rn 2 ~ P2N2, of H2 molecules released to the gas phase. Note 
that according to this formulation F — nnvna^, p\ — 1/tu 
and a = N~ 2 vf{T, a, SE)-yl. 

We will now consider the steady state conditions, where 
N\ = N 2 = 0. In this case iVi can be extracted from Eq. 
@, giving rise to Ni = [-pi + {p'{ + 8aF) 1/2 ]/(4a), where 
the unphysical negative solution is discarded. The steady 
state condition and Eq. (3b) imply Rn 2 = P2N2 = aN 2 . By 
substituting the expression for Ni into this equality we find 
an exact formula for the H2 production rate of a single grain 



pj-pi(pi + SaF) 1/2 +4aF 
8a 



(4) 



Note that the H2 desorption rate p 2 , does not affect the 
steady state production rate Rn 2 ■ However, it will affect the 
number of H2 molecules on thegrain N 2 in the steady state. 
We now evaluate expression (|4j) in two limiting cases. The 
first case is when the adatom desorption rate is negligible 
compared to their recombination rate on the surface. This 
is the limit of pi 2 <C aF. We can neglect the first two terms 
in the numerator on the right hand side of Eq. (ft), finding 



Rn 2 = ^F; 



pi -C aF, 



(•>) 



namely, all H atoms that attach to the surface recombine 
and desorb as H2 molecules. In the other limit, pi 2 aF, 
adatom desorption is important and we expect the produc- 
tion rate Rn 2 to be dependent on p\ . Indeed, expanding the 
square root in Eq. (0) according to (1+a;) 1 ^ 2 ~ l+a;/2 — x 2 /8 
with x — (8aF)/pi , we obtain: 



— F 2 - 

pi 1 



pi > aF. 



(6) 



The important result of Eqs. (g) and (J6J) , which is the focus 
of this note, is that the production rate Rn 2 can be either 
linear or quadratic in the adsorption rate F of H atoms, de- 
pending on the conditions stated above. Eq. (p|) also predicts 
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the dependence of Ru 2 on the recombination rate a and the 
adatom desorption rate pi when the latter is significant. 

To make contact with Eqs. ([j]) and (Q) we substitute 
F = n H t>HO"£, Pi — 1/tfj and a = N~ 2 vf(T,a,5E)jl, and 
multiply the result by rie, the grain density. Thus, Eq. (^|) 
is transformed into Eq. (fil) apart from a factor r\ which is in 
fact e qual t o unity at the conditions studied by Hollenbach 
et al. (1971). Similarly, Eq. (^) is identical to Eq. (H), apart 
from a factor of 1/2 that arises due to a different definition 
of a. 



3 DISCUSSION 

It has been shown analytically that the expressions to calcu- 
late the H2 formation rate given in Eqs. (0) and (^) are cor- 
rect in two different regimes. Eq. (|l|), introduced by Hollen- 
bach et al. (1971) almost thirty years ago, holds when high 
flux, or high mobility of H adatoms, makes it very unlikely 
that H adatoms evaporate before they recombine. Specifi- 
cally, it is the appropriate expression to use when aF S> pi 2 - 
The ass umption of a high mobi li ty has been experiment ally 
proven (Pirronello et al. 1997a; Pirronello et al. 1997b) to 
be inappropriate for the case in which a more realistic grain 
analogue surface is used rather than a monocrystalline one. 
Experimental data suggest that H adatoms have limited mo- 
bility and, in the low coverage regime, H2 is not readily 
formed at the lowest temperature. Such a result had been 



a denser stage, the conditions will probably tend to favour 
Eq. ([l]) due to the increase of the gas phase density, and to 
the decrease of the grain temperature. In passing, we note 
that one must be aware that the processes described above 
of H2 formation on bare refractory grains are influenced by 
at least two other processes, i.e. the decrease of H atom den- 
sity in the gas phase (due to the ongoing conversion into H2) 
and the accretion of an icy mantle. In the latter case, the 
formation rate of H2 might be quite different than the one 
on a sparsely covered silicate surface due to the different 
depth distributions of adsorption sites and hence the differ- 
ent values of sticking and mobility of H adatoms. It would 
be very useful to have laboratory measurements of H2 for- 
mation under these conditions. 

Finally, when a full dense cloud environment is devel- 
oped, Eq. (bp should apply again, because of the low density 
of H in the gas phase and because of the competition be- 
tween H2 and H adatoms in occupying available adsorption 
sites. Moreover, under these conditions, the grain surface 
may be covered by ice where binding energies for H2 are 
slightly larger than for H. In this case, H adatoms will likely 
be far apart from each other and will have to undergo a 
long migration before encountering each other, which in turn 
favours Eq. (^). If we consider evolution from a dense toward 
a diffuse stage the applicability of Eqs. (|l|) and (^) will of 
course take place in the reverse order. Quantitative stud- 
ies of time dependent chemical models of interstellar clouds 



have been done under both static conditions ( 1'Hendecourt 



already | obtained theoretically by Smoluchowski ( 1981 ). In Allamandola fc Greenbergl985|; [Hasegawa fc Herbst 1993| 



Bergin, Langer fc Goldsmith 1995[), a nd dynamically evolv- 
ing conditions ( Ch£xn|ej^eJ^al^T988a ; 5hamle^^d^n)88bj 



his calculations of H2 formation on amorphous water, he ob- 
tained that the onset temperature for mobility and recom- 
binatio n was about 18K, a value significantly higher than 10 



Ncjad, Williams fc Charnlcy 1990]; Prasad, Heere & Tarafdar 



1991| ; [Rawlings et al. 1992| ; jghalabioa fc Groeuberg 1995^ 



K meas ured by Pirronello et al. for olivine. In dense, cold 
clouds, Smoluchowski's value would be too high and would 
make H recombination too infrequent. The express ion given 
in Eq. (^|) introduced by Pirronello et al. (1997b), can be 
applied in the opposite limit or whenever the physical con- 
ditions maintain a low flux or low mobility of H adatoms, 
with respect to their desorption rate. 



IMPLICATIONS FOR THE INTERSTELLAR 
MEDIUM 

The results on the rate o f formation of molecula r hydrogen 



obtained ex perimentally (Pirronello et al. 1997a; Pirronello 
et al. 1997b), and also shown to hold on theoretical grounds 
in this communication, can be used to speculate on how this 
fundamental process might proceed during the dynamical 
evolution of an interstellar cloud. 

As is well known, a typical interstellar cloud crosses dif- 
fuse and dense stages according to its energy balance with 
the surroundings. When a cloud loses more energy than it 
gains from the environment, it contracts; in the opposite 
limit, it expands. In the evaluation of the rate of formation 
of molecular hydrogen, the criteria of applicability of Eqs. 
(jjj) and (Q) will then depend on the particular stage the 
interstellar cloud is going through. In diffuse clouds condi- 
tions favour Eq. (^|). This is because of the low flux of H 
atoms (due to the low gas phase density) and of the rela- 
tively high grain temperature (hence low H residence time 
on the grain surface) . If such a diffuse cloud evolves towards 



The qualitative framework we have just given has to be con- 
firmed by quantitatively incorporating it into calculations 
describing the evolution of an interstellar cloud. 
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